Show all our lenses

Dependencies: ISET3d-v4, ISETCam, isetlens
D. Cardinal, 2022
Note: As of March, 2022, there is a bug in the 2021 version of Matlab, so that the images don't properly render in the LiveScript. So you can either: run this as a Script (using Right-Click from the file tree), or snag an early release of 2022a, or simply look at the HTML version that we've exported from a properly running system.

Initialize ISET and Docker

ieInit;
if ~piDockerExists, piDockerConfig; end

Choose Options

sensorSize =24; % (diagonal in mm)
useScene = 'chessSet'; % one of our default scenes
raysPerPixel = 128;
sensorResolution = 400;
fprintf('Lens Inventory using Scene: %s\n', useScene);
Lens Inventory using Scene: chessSet
fprintf('Sensor Size: %d, rayPP %d, resolution %d\n',...
sensorSize,raysPerPixel,sensorResolution);
Sensor Size: 24, rayPP 128, resolution 400

Read the selected Scene and add a Light

thisR = piRecipeDefault('scene name',useScene);
Read 9 materials. Read 3 textures. ***Scene parsed.
lightName = 'from camera';
ourLight = piLightCreate(lightName,...
'type','distant',...
'cameracoordinate', true);
recipeSet(thisR,'lights', ourLight,'add');

Set render quality

% Set resolution for speed or quality.
thisR.set('film resolution',[sensorResolution sensorResolution]); % 2 is high res. 0.25 for speed
thisR.set('rays per pixel',raysPerPixel); % 128 for high quality
% Pick out a bit of the image to look at. Middle dimension is up.
% Third dimension is z. For the chess set we use a from/to that put the ruler in the
% middle. The in focus is about the pawn or rook.
switch(useScene)
case 'chessSet'
thisR.set('from',[0 0.14 -0.7]); % Get higher and back away than default
thisR.set('to', [0.05 -0.07 0.5]); % Look down default compared to default
case 'cornell_box'
%
end
%thisR.set('aperture diameter',1); % thisR.summarize('all');

Determine the range of object depths in our scene

Note: Getting depth will use pinhole optics, but we use an omni camera with lenses for the rest.
depthRange = thisR.get('depth range');
Complete render took: 3.94e+01 seconds.*** Rendering time for ChessSet: 39.4 sec *** exr2bin done.
%depthRange = [0.1674, 3.3153]; % Chess set distances in meters

Create cameras that use each lens

We get a list of all our lenses, then iterate over it creating a camera object using that lens, and render the same scene through each of them.
lensFiles = lensC.list('quiet', true);
images = []; % our output images
for ii = 1:numel(lensFiles)
lensFileName = lensFiles(ii).name;
lensFile = fullfile(lensFiles(ii).folder, lensFileName);
fprintf('Using lens #%d: %s\n',ii, lensFile);
thisR.camera = piCameraCreate('omni','lensFile',lensFile);
thisR.set('focal distance',mean(depthRange));
% depending on the film size we'll get
% different fields of view
thisR.set('film diagonal',sensorSize);

Render and display

piWrite(thisR);
try
oi = piRender(thisR,'render type',{'radiance', 'depth'});
oi = oiSet(oi,'name',sprintf('Lens: %s ',lensFileName));
fprintf('----------------------------\n');
fprintf('Rendered using: %s',lensFileName);
img = oiShowImage(oi,-1);
images{ii} = img; % for creating a montage
f = figure('Name',lensFileName,'MenuBar','none','ToolBar','none');
imshow(img,'Border','tight');
drawnow; %doesn't seem to help
fprintf('----------------------------\n');
catch
warning('Lens: %s failed to render.\n', lensFile);
end
end
Using lens #1: B:\iset\isetlens\data\lens\2ElLens.json
Complete render took: 1.32e+01 seconds.*** Rendering time for ChessSet: 13.2 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\2ElLens.json
----------------------------
Rendered using: 2ElLens.json
----------------------------
Using lens #2: B:\iset\isetlens\data\lens\2el.XXdeg.100mm.json
Complete render took: 1.34e+01 seconds.*** Rendering time for ChessSet: 13.4 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: 2el.XXdeg.100mm.json
----------------------------
Using lens #3: B:\iset\isetlens\data\lens\2el.XXdeg.12mm.json
Complete render took: 1.18e+01 seconds.*** Rendering time for ChessSet: 11.8 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: 2el.XXdeg.12mm.json
----------------------------
Using lens #4: B:\iset\isetlens\data\lens\2el.XXdeg.50mm.json
Complete render took: 1.31e+01 seconds.*** Rendering time for ChessSet: 13.1 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: 2el.XXdeg.50mm.json
----------------------------
Using lens #5: B:\iset\isetlens\data\lens\2el.XXdeg.6.0mm.json
Complete render took: 1.11e+01 seconds.*** Rendering time for ChessSet: 11.1 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: 2el.XXdeg.6.0mm.json
----------------------------
Using lens #6: B:\iset\isetlens\data\lens\dgauss.22deg.100.0mm.json
Complete render took: 1.47e+01 seconds.*** Rendering time for ChessSet: 14.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.100.0mm.json
----------------------------
Using lens #7: B:\iset\isetlens\data\lens\dgauss.22deg.12.5mm.json
Complete render took: 1.06e+01 seconds.*** Rendering time for ChessSet: 10.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.12.5mm.json
----------------------------
Using lens #8: B:\iset\isetlens\data\lens\dgauss.22deg.3.0mm.json
Complete render took: 9.56e+00 seconds.*** Rendering time for ChessSet: 9.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.3.0mm.json
----------------------------
Using lens #9: B:\iset\isetlens\data\lens\dgauss.22deg.50.0mm.json
Complete render took: 1.37e+01 seconds.*** Rendering time for ChessSet: 13.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.50.0mm.json
----------------------------
Using lens #10: B:\iset\isetlens\data\lens\dgauss.22deg.6.0mm.json
Complete render took: 9.72e+00 seconds.*** Rendering time for ChessSet: 9.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.6.0mm.json
----------------------------
Using lens #11: B:\iset\isetlens\data\lens\dgauss.22deg.6.0mm_v3.json
Complete render took: 9.86e+00 seconds.*** Rendering time for ChessSet: 9.9 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: dgauss.22deg.6.0mm_v3.json
----------------------------
Using lens #12: B:\iset\isetlens\data\lens\dgauss.77deg.3.5201mm.json
Complete render took: 9.70e+00 seconds.*** Rendering time for ChessSet: 9.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\isetlens\data\lens\dgauss.77deg.3.5201mm.json
----------------------------
Rendered using: dgauss.77deg.3.5201mm.json
----------------------------
Using lens #13: B:\iset\isetlens\data\lens\endoscope.json
Complete render took: 3.88e+00 seconds.
Warning: Docker did not run correctly
Status: 1 Result: pbrt version 4 (built Mar 4 2022 at 17:58:41) Copyright (c)1998-2021 Matt Pharr, Wenzel Jakob, and Greg Humphreys. The source code to pbrt (but *not* the book contents) is covered by the Apache 2.0 License. See the file LICENSE.txt for the conditions of the license. [1m[31mWarning[0m: Specified aperture radius 0.0025000002 is greater than maximum possible 0.000081041755. Clamping it. [1m[31mError[0m: Unable to trace ray from film to scene for thick lens approximation. Is aperture stop extremely small? *** Rendering time for ChessSet: 3.9 sec ***
Warning: Lens: B:\iset\isetlens\data\lens\endoscope.json failed to render.
Using lens #14: B:\iset\isetlens\data\lens\endoscopeFixed.json
Complete render took: 3.90e+00 seconds.
Warning: Docker did not run correctly
Status: 1 Result: pbrt version 4 (built Mar 4 2022 at 17:58:41) Copyright (c)1998-2021 Matt Pharr, Wenzel Jakob, and Greg Humphreys. The source code to pbrt (but *not* the book contents) is covered by the Apache 2.0 License. See the file LICENSE.txt for the conditions of the license. [1m[31mWarning[0m: Specified aperture radius 0.0025000002 is greater than maximum possible 0.0020584606. Clamping it. [1m[31mError[0m: Unable to trace ray from film to scene for thick lens approximation. Is aperture stop extremely small? *** Rendering time for ChessSet: 3.9 sec ***
Warning: Lens: B:\iset\isetlens\data\lens\endoscopeFixed.json failed to render.
Using lens #15: B:\iset\isetlens\data\lens\fisheye.87deg.100.0mm.json
Complete render took: 1.33e+01 seconds.*** Rendering time for ChessSet: 13.3 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\fisheye.87deg.100.0mm.json
----------------------------
Rendered using: fisheye.87deg.100.0mm.json
----------------------------
Using lens #16: B:\iset\isetlens\data\lens\fisheye.87deg.12.5mm.json
Complete render took: 1.23e+01 seconds.*** Rendering time for ChessSet: 12.3 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\fisheye.87deg.12.5mm.json
----------------------------
Rendered using: fisheye.87deg.12.5mm.json
----------------------------
Using lens #17: B:\iset\isetlens\data\lens\fisheye.87deg.3.0mm.json
Complete render took: 1.02e+01 seconds.*** Rendering time for ChessSet: 10.2 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\fisheye.87deg.3.0mm.json
----------------------------
Rendered using: fisheye.87deg.3.0mm.json
----------------------------
Using lens #18: B:\iset\isetlens\data\lens\fisheye.87deg.50.0mm.json
Complete render took: 1.40e+01 seconds.*** Rendering time for ChessSet: 14.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\fisheye.87deg.50.0mm.json
----------------------------
Rendered using: fisheye.87deg.50.0mm.json
----------------------------
Using lens #19: B:\iset\isetlens\data\lens\fisheye.87deg.6.0mm.json
Complete render took: 1.16e+01 seconds.*** Rendering time for ChessSet: 11.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\fisheye.87deg.6.0mm.json
----------------------------
Rendered using: fisheye.87deg.6.0mm.json
----------------------------
Using lens #20: B:\iset\isetlens\data\lens\gullstrand.json
Complete render took: 1.14e+01 seconds.*** Rendering time for ChessSet: 11.4 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\gullstrand.json
----------------------------
Rendered using: gullstrand.json
----------------------------
Using lens #21: B:\iset\isetlens\data\lens\microlens.json
Complete render took: 3.77e+00 seconds.
Warning: Docker did not run correctly
Status: 1 Result: pbrt version 4 (built Mar 4 2022 at 17:58:41) Copyright (c)1998-2021 Matt Pharr, Wenzel Jakob, and Greg Humphreys. The source code to pbrt (but *not* the book contents) is covered by the Apache 2.0 License. See the file LICENSE.txt for the conditions of the license. [1m[31mWarning[0m: Specified aperture radius 0.0025000002 is greater than maximum possible 0.0000010000001. Clamping it. [1m[31mError[0m: Unable to trace ray from scene to film for thick lens approximation. Is aperture stop extremely small? *** Rendering time for ChessSet: 3.8 sec ***
Warning: Lens: B:\iset\isetlens\data\lens\microlens.json failed to render.
Using lens #22: B:\iset\isetlens\data\lens\petzval.12deg.100.0mm.json
Complete render took: 1.30e+01 seconds.*** Rendering time for ChessSet: 13.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\petzval.12deg.100.0mm.json
----------------------------
Rendered using: petzval.12deg.100.0mm.json
----------------------------
Using lens #23: B:\iset\isetlens\data\lens\petzval.12deg.12.5mm.json
Complete render took: 9.71e+00 seconds.*** Rendering time for ChessSet: 9.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\petzval.12deg.12.5mm.json
----------------------------
Rendered using: petzval.12deg.12.5mm.json
----------------------------
Using lens #24: B:\iset\isetlens\data\lens\petzval.12deg.3.0mm.json
Complete render took: 9.20e+00 seconds.*** Rendering time for ChessSet: 9.2 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\petzval.12deg.3.0mm.json
----------------------------
Rendered using: petzval.12deg.3.0mm.json
----------------------------
Using lens #25: B:\iset\isetlens\data\lens\petzval.12deg.50.0mm.json
Complete render took: 1.27e+01 seconds.*** Rendering time for ChessSet: 12.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\petzval.12deg.50.0mm.json
----------------------------
Rendered using: petzval.12deg.50.0mm.json
----------------------------
Using lens #26: B:\iset\isetlens\data\lens\petzval.12deg.6.0mm.json
Complete render took: 9.65e+00 seconds.*** Rendering time for ChessSet: 9.7 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\petzval.12deg.6.0mm.json
----------------------------
Rendered using: petzval.12deg.6.0mm.json
----------------------------
Using lens #27: B:\iset\isetlens\data\lens\reversed.telephoto.42deg.27mm.json
Complete render took: 1.32e+01 seconds.*** Rendering time for ChessSet: 13.2 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\isetlens\data\lens\reversed.telephoto.42deg.27mm.json
----------------------------
Rendered using: reversed.telephoto.42deg.27mm.json
----------------------------
Using lens #28: B:\iset\isetlens\data\lens\reversed.telephoto.42deg.3.5201mm.json
Complete render took: 1.00e+01 seconds.*** Rendering time for ChessSet: 10.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\isetlens\data\lens\reversed.telephoto.42deg.3.5201mm.json
----------------------------
Rendered using: reversed.telephoto.42deg.3.5201mm.json
----------------------------
Using lens #29: B:\iset\isetlens\data\lens\reversed.telephoto.77deg.3.5201mm.json
Complete render took: 9.96e+00 seconds.*** Rendering time for ChessSet: 10.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\isetlens\data\lens\reversed.telephoto.77deg.3.5201mm.json
----------------------------
Rendered using: reversed.telephoto.77deg.3.5201mm.json
----------------------------
Using lens #30: B:\iset\isetlens\data\lens\telephoto.250mm.json
Complete render took: 3.88e+00 seconds.
Warning: Docker did not run correctly
Status: 1 Result: pbrt version 4 (built Mar 4 2022 at 17:58:41) Copyright (c)1998-2021 Matt Pharr, Wenzel Jakob, and Greg Humphreys. The source code to pbrt (but *not* the book contents) is covered by the Apache 2.0 License. See the file LICENSE.txt for the conditions of the license. [1m[31mError[0m: Coefficient must be positive. It looks focusDistance 1 is too short for a given lenses configuration *** Rendering time for ChessSet: 3.9 sec ***
Warning: Lens: B:\iset\isetlens\data\lens\telephoto.250mm.json failed to render.
Using lens #31: B:\iset\isetlens\data\lens\tessar.22deg.100.0mm.json
Complete render took: 1.45e+01 seconds.*** Rendering time for ChessSet: 14.5 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: tessar.22deg.100.0mm.json
----------------------------
Using lens #32: B:\iset\isetlens\data\lens\tessar.22deg.12.5mm.json
Complete render took: 1.10e+01 seconds.*** Rendering time for ChessSet: 11.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: tessar.22deg.12.5mm.json
----------------------------
Using lens #33: B:\iset\isetlens\data\lens\tessar.22deg.3.0mm.json
Complete render took: 9.60e+00 seconds.*** Rendering time for ChessSet: 9.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: tessar.22deg.3.0mm.json
----------------------------
Using lens #34: B:\iset\isetlens\data\lens\tessar.22deg.50.0mm.json
Complete render took: 1.36e+01 seconds.*** Rendering time for ChessSet: 13.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: tessar.22deg.50.0mm.json
----------------------------
Using lens #35: B:\iset\isetlens\data\lens\tessar.22deg.6.0mm.json
Complete render took: 9.96e+00 seconds.*** Rendering time for ChessSet: 10.0 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: tessar.22deg.6.0mm.json
----------------------------
Using lens #36: B:\iset\isetlens\data\lens\wide.40deg.100.0mm.json
Complete render took: 1.11e+01 seconds.*** Rendering time for ChessSet: 11.1 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.40deg.100.0mm.json
----------------------------
Using lens #37: B:\iset\isetlens\data\lens\wide.40deg.12.5mm.json
Complete render took: 1.13e+01 seconds.*** Rendering time for ChessSet: 11.3 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.40deg.12.5mm.json
----------------------------
Using lens #38: B:\iset\isetlens\data\lens\wide.40deg.3.0mm.json
Complete render took: 9.06e+00 seconds.*** Rendering time for ChessSet: 9.1 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.40deg.3.0mm.json
----------------------------
Using lens #39: B:\iset\isetlens\data\lens\wide.40deg.50.0mm.json
Complete render took: 1.16e+01 seconds.*** Rendering time for ChessSet: 11.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.40deg.50.0mm.json
----------------------------
Using lens #40: B:\iset\isetlens\data\lens\wide.40deg.6.0mm.json
Complete render took: 9.52e+00 seconds.*** Rendering time for ChessSet: 9.5 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.40deg.6.0mm.json
----------------------------
Using lens #41: B:\iset\isetlens\data\lens\wide.56deg.100.0mm.json
Complete render took: 1.23e+01 seconds.*** Rendering time for ChessSet: 12.3 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.100.0mm.json
----------------------------
Using lens #42: B:\iset\isetlens\data\lens\wide.56deg.12.5mm.json
Complete render took: 1.12e+01 seconds.*** Rendering time for ChessSet: 11.2 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.12.5mm.json
----------------------------
Using lens #43: B:\iset\isetlens\data\lens\wide.56deg.3.0mm.json
Complete render took: 9.40e+00 seconds.*** Rendering time for ChessSet: 9.4 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.3.0mm.json
----------------------------
Using lens #44: B:\iset\isetlens\data\lens\wide.56deg.50.0mm.json
Complete render took: 1.25e+01 seconds.*** Rendering time for ChessSet: 12.5 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.50.0mm.json
----------------------------
Using lens #45: B:\iset\isetlens\data\lens\wide.56deg.6.0mm.json
Complete render took: 1.06e+01 seconds.*** Rendering time for ChessSet: 10.6 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.6.0mm.json
----------------------------
Using lens #46: B:\iset\isetlens\data\lens\wide.56deg.6.0mm_v3.json
Complete render took: 1.08e+01 seconds.*** Rendering time for ChessSet: 10.8 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.
----------------------------
Rendered using: wide.56deg.6.0mm_v3.json
----------------------------
Using lens #47: B:\iset\isetlens\data\lens\wide.77deg.4.38mm.json
Complete render took: 9.83e+00 seconds.*** Rendering time for ChessSet: 9.8 sec *** exr2bin done.exr2bin done.exr2bin done.exr2bin done.Decoding lens parameters from lens file name B:\iset\iset3d-v4\data\lens\wide.77deg.4.38mm.json
----------------------------
Rendered using: wide.77deg.4.38mm.json
----------------------------

Show a Captioned Montage for Reference

For filename captions to be legible we need to make sure we have enough room
minTileSize = 400; thumbnailSize = 400;
boxSize = max(minTileSize, sensorResolution);
textPosition = [10 10; boxSize-30 boxSize-30];
for ii = 1:numel(images)
if minTileSize > sensorResolution && ~isempty(images{ii})
imagesR{ii} = imresize(images{ii},[boxSize boxSize]);
else
imagesR{ii} = images{ii};
end
imagesC{ii} = insertText(imagesR{ii}, textPosition, lensFiles(ii).name, ...
'FontSize', 36); %#ok<*SAGROW>
end
montage(imagesC,'ThumbnailSize',[thumbnailSize thumbnailSize]);